Self-interacting Scalar Field Trapped in a Randall-Sundrum Braneworld: The 

Dynamical Systems Perspective 
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We apply the dynamical systems tools to study the linear dynamics of a self-interacting scalar 
field trapped on a Randall-Sundrum brane. The simplest kinds of self-interaction potentials are 
investigated: a) constant potential, and b) exponential potential. It is shown that the dynamics 
of the Randall-Sundrum model significantly differs from the standard four-dimensional behavior 
at early times: in all cases of interest the (singular) empty universe is the past attractor for every 
trajectory in phase space, meanwhile the kinetic energy-dominated solution is always a saddle critical 
point. The late-time dynamics is not affected by the brane effects. 
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I. INTRODUCTION 

Recent observations from the Wilkinson Microwave 
Anisotropy Probe (WMAP) [H offer strong supporting 
evidence in favor of the inflationary paradigm. In the 
most simple models of this kind, the energy density of the 
universe is dominated by the potential energy of a single 
(inflaton) scalar field that slowly rolls down in its self- 
interaction potential Restrictions imposed upon the 
class of potentials which can lead to realistic inflationary 
scenarios, are dictated by the slow-roll approximation, 
and hence, the result is that only sufficiently flat poten- 
tials can drive inflation. In order for the potential to be 
sufficiently flat, these conventional inflationary models 
should be fine-tuned. This simple picture of the early- 
time cosmic evolution can be drastically changed if one 
considers models of inflation inspired in "Uniflcd Theo- 
ries" like the Super String or M-theory. One of the most 
appealing models of this kind is the Randall-Sundrum 
braneworld model of type 2 (RS2) [3]. In this model a 
single codimension 1 brane with positive tension is em- 
bedded in a five-dimensional anti-de Sitter (AdS) space- 
time, which is infinite in the direction perpendicular to 
the brane. In general, the standard model particles are 
conflned to the brane, meanwhile gravitation can propa- 
gate in the bulk. In the low-energy limit, due to the cur- 
vature of the bulk, the graviton is confined to the brane, 
and standard (four-dimensional) general relativity laws 
are recovered. 

RS2 braneworld models have an appreciable impact on 
early universe cosmology, in particular, for the inflation- 
ary paradigm. In effect, a distinctive feature of cosmol- 



ogy with a scalar fleld confined to a RS2 brane is that 
the expansion rate of the universe differs at high energy 
from that predicted by standard general relativity. This 
is due to a term - quadratic in the energy density - that 
produces enhancing of the friction acting on the scalar 
field. This means that, in RS2 braneworld cosmology, 
inflation is possible for a wider class of potentials than in 
standard cosmology [4] . Even potentials that are not suf- 
ficiently flat from the point of view of the conventional 
inflationary paradigm can produce successful inflation. 
At sufficiently low energies (much less than the brane 
tension), the standard cosmic behavior is recovered prior 
to primordial nucleosynthesis scale (T ~ 1 MeV) and a 
natural exit from inflation ensues as the fleld accelerates 
down its potential 0].^ 

Another interesting feature of this scenario is that the 
inflaton does not necessarily need to decay; it may sur- 
vive through the present epoch in the cosmic evolution. 
Therefore, it may also play the role of the quintessence 
field, which is a necessary ingredient to explain the cur- 
rent acceleration in the expansion of the universe. Such a 
unified theoretical framework for the description of both 
inflaton and quintessence with the help of just one single 
scalar field has been the target of some works (see for 
instance Refs. [1, [13, HH, 111 ) • 

Aim of this letter is to study the dynamics of a such 
scenario - a self-interacting scalar field trapped on a RS2 
brane - by invoking the dynamical systems tools. These 
have been proved useful to retrieve significant informa- 
tion about the evolution of a huge class of cosmological 
models (see for instance the book J^] ) . The simplest self- 
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In this scenario, reheating arises naturally even for potentials 
without a global minimum and radiation is created through grav- 
itational particle production Q and/or through curvaton reheat- 
ing Q. This last ingredient improves the brane "steep" infla- 
tionary picture Q. Other mechanisms such as preheating, for 
instance, have also been explored [^. 
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interaction potentials: a) constant potential, and b) ex- 
ponential potential, are investigated. The latter one rep- 
resents a common functional form for self-interaction po- 
tentials that can be found in higher-order or higher- 
dimensional theories [1^]. These can also arise due to 
non-perturbative effects [l6l |. In addition to the scalar 
field we also consider a background fluid trapped on the 
RS brane. 

In a sense, the present research is a generalization of 
the one reported in [l3|, to include higher-dimensional 
behaviour (in the present case dictated by the RS dy- 
namics). In consequence, for the exponential potential, 
our results will include the ones reported in as a 
particular case. This is true even if the early-time (high- 
energy) dynamics does actually differ from the standard 
(four-dimensional) one, regarding the stability properties 
of the critical points. This inve stig ation is, also, comple- 
mentary to the one in reference fl8l|, where the target was 
the study of the dynamics of a single scalar field trapped 
in a Dvali-Gabadadze-Porrati (DGP) brane. In the latter 
case, however, only the late-time dynamics is modified, 
while the dynamics of early times is not affected in any 
essential way. 

Through the paper we use natural units (SttG = 
Stt/top; — c~ 1). 

II. PHASE SPACE 



This means that the 4D (low-energy) limit of the Randall- 
Sundrum cosmological equations - corresponding to the 
formal limit A — s- cx) - is associated with the value z = 1. 
The high-energy limit A — > 0, on the contrary, corre- 
sponds to z — > 0. However, wc shall exclude the points 
with z = (and their neighborhood) since, in general, 
at 2; = the equations () are undcfeined. The critical 
point associated with z — has to be analyzed sepa- 
ratelly. In any case the physical meaning of this point in 
phase space has to be taken with caution due to the high 
energies associated with it. 

It arises the following constraint: 



Since p„i > 0, then < x'^+y'^ < z. We will be focused on 
expanding FRW universes, so that y > 0. The resulting 
phase space for the RS model is the following: 

* = {ix,y,z) -.0 <x^ +y^ < z, 

-1 < a; < 1,0 < 2/,0 < z < 1}. (7) 

After the above choice of variables of the phase space, 
the system of cosmological equations can be trans- 
lated into the following autonomous system of ordinary 
differential equations: 



The starting point will be the cosmological equations 
for the RS2 model which, for a Friedmann-Robertson- 
Walker (FRW) universe with flat spatial sections, can be 
written in the following way: 
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where pt = P<t, + Pm {P4> = <^^/2 -|- T^, V - scalar field self- 
interaction potential) is the total energy content of the 
brane, A is the brane tension, 7 is the barotropic index of 
the background fluid, and the dot accounts for derivative 
in respect to the cosmic time t. Our aim is to write the 
latter system of second-order (partial) differential equa- 
tions, as an autonomous system of (first order) ordinary 
differential equations. For this purpose we introduce the 
following phase variables: 
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After this choice of variables one can realize that 
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y[2x^ + ^{z~x^-y^)], (9) 



z' = 3(1 - z)[2x2 -I- 7(z - - y% 



(10) 



where the prime denotes derivative with respect to the 
new time variable r = In a. 

It will be helpful to have the parameters of observa- 
tional importance fi^ — p^/iH"^ - the scalar field dimen- 
sionless energy density parameter, and the equation of 
state (EOS) parameter uj^ = (t/f^ - 2V)/{(j) + 2V), writ- 
ten in terms of the variables of phase space: 



x'^+y^, UJ4, 



2 ? 
X — y 

+ y^ 



(11) 



Additionally, the decceleration parameter q = —(I 



'^ + l(^]i^-' + ii^~-'~y')]- (12) 
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TABLE I: Properties of the critical points for the autonomous system (|13p . 
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TABLE IL Eigenvalues for the critical points in table U 
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III. THE CRITICAL POINTS 

As in [ill we study here two concrete yet generic cases. 
First we consider a constant potential V — Vq, and then 
we study the exponential potential V ~ Vo exp (— at/)). 
For more general classes of potentials we have to rely 
on a quite different approach (see, for instance reference 
[1^). This will be the subject of future research. 

A. The Constant Potential V = Vo 

The autonomous system of ordinary differential equa- 
tions ((8] fT0|l simplifies to: 

x' = -3x + ^ {^~~~^ ^[^^^ + ~ y^)], 

z' = i{l - z)[2x^ + -i{z - - y^)]. (13) 

The critical points of the autonomous system of equa- 
tions (fT5)) are summarized in table [H while table HIl shows 
the eigenvalues of the corresponding linearization (Jaco- 
bian) matrices. 

Critical points of the autonomous system (|13p can be 
divided into two sets according to whether z ^ 1 or 
z = 1. The latter set of points is associated with four- 
dimensional behavior, meanwhile for the critical points 
with z ^ 1 five-dimensional effects are present. There is 
just one critical point with z ^ \ (point P4 in table [l|: 




FIG. 1: Trajectories in phase space for given initial conditions 
for the constant potential V = Vq (upper part of the figure: 
view from top to botton in the left panel, and from botton 
to top in the right panel). The fiux in time r is shown in the 
lower part of the figure. At late times trajectories in phase 
space approach to the plane {x, y, 1) that is associated with 
standard 4D behavior. The late-time (future) attractor is the 
infiationary de Sitter-FRW universe. 



(x, y, z) — (0, z). Strictly speaking it is not, in fact, a 
critical point but a subset of critical points. In this case, 
since z = y"^ (x = 0), then pm = 0, while pT = V, so that 
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the Friedmann equation takes the form: 



3H' 



V I 



V 
2A 



(14) 



This represents the slow-roll Friedmann equation relating 
the Hubble expansion parameter with the potential of 
the inflaton field, modified by the presence of the RS 
brane. For potentials much larger than the brane tension 
(properly the limit when the effects coming from the RS 
brane are appreciable): V ^ \ ^ Hrs = V/V6X, so 
that the early-time (high-energy) expansion rate in the 
Randall-Sundrum model Hjig gets enhanced with respect 
to the general relativity rate Hgr- 



H 



RS 



Hgr. 



(15) 



This is the way brane effects fuel early inflation in the 
RS model. As seen from table HIl one of the eigenvalues 
of the linearization matrix for P4 vanishes, so that it is a 
non- hyperbolic point (set of points), and the linear anal- 
ysis is not conclusive in this case. However the analysis 
of the trajectories in phase space reveals the true struc- 
ture of the phase space. In the figure [T] it is shown how 
the trajectories in phase space emerge from the point 
S = {x,y,z) = (0,0,0) - the empty (Misner-RS) uni- 
verse - meaning that this is the past attractor of the RS 
cosmological model. This can be verified by expanding 
the equations ((8l fT0|) in the neighborhood of 5 = (0, 0, 0), 
i. e., by replacing y + e, z— >0 + e, while keeping 
X — (e - a very small perturbation). The result is that 
the perturbation e increases exponentially with time r: 
e(T) = eo exp (3jt), so that this point is unstable. 

A cautionary note is neccessary: we want to notice 
that the points with z = have been removed from the 
phase space 4* since, in general, at z = the autonomous 
system of equations (|8l lT0)) blows up due to our choice of 
phase space variables. For that reason the point S does 
not appear in table HI 

There are three additional isolated critical points cor- 
responding to the autonomous system (|13p , which are as- 
sociated with the standard 4D limit of the theory z — 1. 
The first of these critical points is the matter-dominated 
solution: 3iJ^ = p„i (point Pi = (0,0,1) in table |T|. As 
seen from table [H] it is always a saddle point in phase 
space. The second one: — (±1,0,0), is the kinetic 
energy-dominated solution {3H^ — 0^/2), which is also 
a saddle point in phase space. At first sight this is an 
unexpected result since, in standard general relativity, 
the kinetic energy-dominated solution is always a source 
(past attractor) in phase space. However, a closer inspec- 
tion of the RS model reveals that this point belongs in 
the high-energy phase, which is the one that is modified 
by the brane contribution, hence it is not such an unex- 
pected result: five-dimensional contributions modify the 
structure of the phase space at high-energies. 






FIG. 2: Phase trajectories for given initial data, for the expo- 
nential potential V — Vo exp {—acj)) (upper part of the figure: 
view from top to botton in the left panel, and from botton 
to top in the right panel). The corresponding flux in time r, 
is shown in the panels in the lower part of the figure. The 
exponent has been chosen to be a = 2 so that the scalar 
field-dominated solution is the late-time attractor. 



The late-time attractor is the inflationary de Sitter- 
FRW solution = V, which corresponds to the point 
P3 = (0, 1, 1) in table H] (see also the figure [T]). 



B. The exponential Potential V — Vo exp {—a(j>) 

The autonomous system of ordinary differential equa- 
tions ((8] fT0|) takes the following form: 



/3 2 o 
X' = \ / -ay — 3x + 



3 /2- z 



x[2x'^ + 7(z - - y"^)], 



3 /2 



y = \l-axy + - 



y[2x'^ + "f{z — x^ ~ y^ 



z' = 3(1 - z)[2x^ + 7(z -x^ - y^)]. 



(16) 



The critical points of (fT6|) are summarized in table IIIll 
while the corresponding eigenvalues are shown in table 
IIVI In this case there is only one critical point that can 
be associated with higher-dimensional behavior (z 7^ 1): 
S — {x,y,z) = (0,0,0), instead of a whole set of points 
as in the former case. Actually, at the critical point x' — 
y' — z' = 0, therefore, from the third equation in (I16p it 
follows that 
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TABLE III: Critical points for the autonomous system (|16|l . 
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1 


All a, 7 
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Pi 


V 2 Q 
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> 37 


37 
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-1 + 


37 
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TABLE IV: Eigenvalues for the critical points in table lllll In the row corresponding to the last critical point P4, for compactness, 
we have introduced the notation II = -^7(2 — j){24{'y /a)^ — 97 + 2). 



p 


X 


V 


z 


Ai 


A2 


A3 



Pi 1 -f(2-7) -37 37/2 

P^ ±1 1 \l\{a - yfl) -6 3(2 - 7) 

^3 



\/6 



1 -f[(2-7) + n] -|[(2-7)-n] -37 




FIG. 3: Trajectories in phase space for given initial data, for 
the exponential potential V = exp (— a(^) (upper part of 
the figure), and the corresponding flux in time r (lower part 
of the figure). The exponent now has been set a = 4 so that 
the matter-scaling solution is the late-time attractor in this 
case. 



z' = ^ + (17) 

meanwhile, simultaneously, 

a:' = ^ X — —p^iP' — 3x, 

y' = - y/^axy = ^ x = y ^ z ^ 0. (18) 

We have to recall, however, that the points with coor- 
dinate z = have been removed from the phase space. 
This is due to technical reasons related with our choice 
of coordinates. For that reason, as in the former case, 
the point S does not appear in table [TTl 

To understand what happens at the point {x,y,z) = 
(0, 0, 0), it is recommended to expand the equations pB]) 
in its neighborhood, i. e., to replace x^O + e, y^O + e. 
As before one gets that the (homogeneous) perturbation 
e grows exponentially in time r: e(T) — eoexp(37T), so 
that this point is a source (past attractor) in phase space. 

There are four additional (isolated) critical points 
which are associated with the 4D limit. A complete study 
of these critical points (see table IIIip. for the standard 4D 
case (z = 1), can be found in the well known reference 
fll\ . However, at early times, the results of [l3| differ a 
little bit from the ones reported here, due to the modifi- 
cations produced by the RS brane effects on the cosmic 
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Asymptotic solution 




/ / / 

i 



Asymptotic solution 




FIG. 4: Projections of the phase trajectories onto the plane 
{x, 3/, 2 = 1) - associated with low-energy (late-time) behavior 
- for the exponential potential V = Vo exp {—a(j>) (a = 2 in 
the upper figure, while a = 4 in the lower figure). Note 
that, although the kinetic energy-dominated/matter-scaling 
solution (upper/lower part of the figure respectively) appears 
to be the past attractor, it is just a "mirage" - an effect of 
the projection - due to the ignorance of the brane effects. 



dynamics. Our main results can be summarized as fol- 
lows (see tables |lll] and IIVI) : 

• Va,7 (0 < 7 < 2); The past attractor in the phase 
space is the empty universe (point S that does 
not appear in table IIIip (figure [2]) . The matter- 
dominated solution (point Pi in tabic , and the 
kinetic energy-dominated solutions (points ^2^)^ 
are always saddle points in phase space. ^ 

• < 87; The scalar field-dominated solution (point 
P4) is the late-time attractor (see figure [2]). 

• 37 < < 6; The scalar field-dominated solution 
P3 is a saddle point. Whenever 24(7/0;)^ is big- 
ger/lower than 97—2, the matter-scaling solution 
(point P4 in table IIIip is a stable node/spiral (see 
figures [3] and SI . 



> 6; The matter-scaling solution P4 is a stable 
spiral (see figures [3] and 0]) . 



IV. RESULTS AND DISCUSSION 

From the analysis in the former sections, the following 
important results can be summarized: 

• For the constant potential V — Vq, early infla- 
tion fuelled by the potential energy of the scalar 
field trapped in the Randall-Sundrum brane (equa- 
tion (dH)), is a (non-isolated) critical point in phase 
space. ^ The late-time attractor is the (inflationary) 
de Sitter space. 

• For the exponential potential V — Vo exp (— at/)), 
early inflation is not even a critical point. This 
means that, for this potential, inflation is not 
generic even if the higher-dimensional effects in 
the Randall-Sundrum scenario favor early inflation. 
Depending on the values of the free parameters a 
and 7, the late-time attractor can be either the 
scalar field-dominated solution (a^ < 87, infla- 



tionary if 
{a^ > 87). 



< 2), or the matter-scaling solution 



• For both potentials the past (high-energy) attrac- 
tor is the (singular) empty universe (2 = 
Pt — 0)- However, given the high energies implied, 
the physical significance of this point is suspicious. 
The kinetic energy-dominated solution is always a 
saddle point in phase space. 

The latter result marks a difference with respect to the 
well known results within general relativity plus a min- 
imally coupled scalar field [l3|, according to which the 
kinetic energy-dominated solution is the past attractor 
for V = Vq, while, for the exponential potential it can 
be either a past attractor, or a saddle point. The ori- 
gin of this difference is in the high energy modificatios of 
the dynamics caused by the brane effects in the Randall- 
Sundrum scenario, leading to the empty universe being 
always the past attractor.'* 

Another interesting result has to do with the fact 
that, for the constant potential V = Vq, the early-time 
inflation driven by the potential energy of the scalar 
field - properly the cosmological constant - represents 
a critical subset in the phase space of the RS model: 
S = {0,y/z, z), 



^ As it has been properly noted before, due to the influence of the 
brane effects, the early-times dynamics is modified with respect 
to the standard 4D general relativity situation. 



^ As aheady noted, it is in fact a critical set or, strictly speaking, 

a set of critical points: S = (0, y'z, 2), 2 E]0, 1]. 
* As already noted, this critical point has to be taken with some 

care due to the very high energies associated with it. One has 

to be sure of the validity of the Randall-Sundrum model at such 

large energies. 
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^ 3i/2^yo(l + ^). (19) 

For sufficiently large values of the brane cosmological con- 
stant A = Vb ^ A, the Hubble expansion rate is enhanced 
with respect to the standard general relativity (GR) rate: 
H = Vb/TeA. 

This subset of critical points could be associated with 
early inflation. The fact that this solution is a subset 
of critical points in phase space, means that, for V — 
Vb, early inflation described by equation (fTO|) . is quite a 
generic property of the cosmological RS model. For the 
exponential potential, on the contrary, early inflation is 
not even a critical point in the phase space, meaning 
that the genericity of inflation is dependent on the kind 
of self-interaction potential under consideration. 

The phase space pictures in the figures [TJ [51 and El re- 
veal the actual behavior of the RS dynamics: trajectories 
in phase space depart from the (singular) empty universe, 
possibly related with the initial (Big-Bang) singularity, 
and, at late times, approach to the plane (x, y, 1), which 
is associated with standard four-dimensional behavior. 



V. CONCLUSION 

The dynamics of RS models departs from the one in 
general relativity at high energies (early times). In par- 
ticular, a dynamical systems approach to the subject re- 
veals that the kinetic energy-dominated solution, which 
is usually the past attractor for every trajectory in the 
phase space, is replaced by the empty universe. The ki- 
netic energy-dominated phase represents always a saddle 

^ See, for instance, the recent reference [2ll |. 



point. At late times gravity is effectivelly confined to the 
brane so that the late-time dynamics is not modified with 
respect to the standard general relativity case. 

For the constant potential case there are two critical 
points - properly a subset of non-isolated critical points, 
plus an isolated critical point - that are linked with in- 
flationary behavior. This entails a good possibility for 
this model to account for an unified description of early 
inflation and of the late-time accelerated expansion of 
the universe, with the help of a single scalar field. This 
possibility, whithin the frame of the Randall- Sundrum 
model, has been studied, in other contexts, for instance, 
in reference (20| . However, for the exponential poten- 
tial the situation is quite different: there are no criti- 
cal points that could be associated with early inflation. 
Hence, the possibility for an unified description of infla- 
tion and quintessence within RS braneworld model is not 
as generic as thought. 

It can be of insterest to investigate the present scenario 
for arbitrary self- interaction potentials, so that we will 
be able to reach to more generic results. This task would 
entail a different approach than the one undertaken in 
this paper, ^ so that we leave it for future work. 
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